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Abstract 

I discuss the role played by multiple partonic interactions (MPI) in 
the early stages of relativistic heavy-ion collisions, for which a weak- 
coupling QCD description is possible. From the Color Glass Conden- 
sate, through the Glasma and into the Quark-Gluon-Plasma phase, MPI 
are at the origin of interesting novel QCD phenomena. 

1 Introduction 

Relativistic heavy-ion collisions involve such large parton densities, that they are reactions where 
multiple partonic interactions (MPI) abound, and in which those can be investigated. Through 
most of the stages of a high-energy heavy-ion collision, MPI are not only important but cru- 
cial, and without their understanding, no robust QCD-based description of the collision can be 
achieved. During the different phases that the system goes through, from the initial nuclear wave 
functions, through the pre-equilibrium state just after the collision, and into the following ther- 
malized quark-gluon plasma (QGP) and hadronic phases, MPI are at the origin of most interesting 
phenomena. 

However, one may wonder what can be described with first-principle weak-coupling QCD 
calculations. It has been proposed that the early stages of the heavy-ion collision should be, per- 
haps until the QGP phase. The saturation of the initial nuclear wave functions, and the multipar- 
ticle production from the decay of strong color fields are phenomena which have been addressed 
by weak-coupling methods, as well as the quenching of hard probes via QGP-induced energy 
loss. In those calculations, MPI are characterized by momentum scales which, if hard enough, 
justify a weak coupling analysis. 

In the Color Glass Condensate (CGC) picture of the nuclear wave function, the saturation 
scale Q s characterizes which quantum fluctuations can be treated incoherently and which cannot; 
in the glasma phase right after the collision of two CGCs, l/Q s sets the time scale for the decay 
of the strong color fields; and in the QGP phase, the plasma saturation momentum characterizes 
what part of the wave function of hard probes is responsible for their energy loss, by becoming 
emitted radiation. In the following, I discuss the role played by MPI in those different stages. 

2 The saturation scale in the nuclear wave function 

The QCD description of hadrons/nuclei in terms of quarks and gluons depends on the process 
under consideration, on what part of the wave function is being probed. Consider a hadron 
moving at nearly the speed of light along the light cone direction x + , with momentum P + . 
Depending on their transverse momentum hp and longitudinal momentum xP + , the partons 
inside the hadron behave differently, reflecting the different regimes of the hadron wave function. 




Fig. 1: Left: diagram in the (fcy = Q 2 , x) plane picturing the hadron/nucleus in the different weakly-coupled regimes. 
The saturation line separates the dilute (leading-twist) regime from the dense (saturation) regime. Right: when scat- 
tering a dilute probe on the hadron/nucleus, both multiple scatterings and saturation of the wave function are equally 
important at small x, when occupation numbers become of order l/a s . 



When probing the (non-perturbative) soft part of the wave function, corresponding to par- 
tons with transverse momenta of the order of Aqcd ~ 200 MeV, the hadron looks like a bound 
state of strongly interacting partons. When probing the hard part of the wave function, corre- 
sponding to partons with /ctS> Aqcd and x< 1, the hadron looks like a dilute system of weakly 
interacting partons. 

The saturation regime of QCD describes the small— x part of the wave function. When 
probing partons that feature fcf> Aqcd, and x<Cl, the effective coupling constant a s log(l/x) 
is large, and the hadron looks like a dense system of weakly interacting partons, mainly gluons 
(called small— x gluons). The larger kf is, the smallest x needs to be to enter the saturation 
regime. As pictured in Fig.l, this means that the separation between the dense and dilute regimes 
is characterized by a momentum scale Q s (x), called the saturation scale, which increases as x 
decreases. 

A simple way to estimate the saturation scale is to equate the gluon-recombination cross- 
section a rec ~ a s /kj, with 1/pr ~ nR 2 /(xf(x,kj,)), the inverse gluon density per unit of 
transverse area. Indeed, when a rec pT ~ 1, one expects recombination not to be negligible 
anymore. This gives: 

^2 _ a s xf(x,Q 2 s ) 

Note that a s (Q 2 ) decreases as x decreases, so for small enough x, one deals with a weakly- 
coupled regime, even though non-linear effects are important. The scattering of dilute partons 
(with kx S> Q s (x)) is described in the leading-twist approximation in which they scatter incoher- 
ently. By contrast, when the partem density is large (kr ~ Q s (x)), partons scatter collectively. 

The Color Glass Condensate (CGC) is an effective theory of QCD [1] which aims at de- 
scribing this part of the wave function. Rather than using a standard Fock-state decomposition, 
it is more efficient to describe it with collective degrees of freedom, more adapted to account for 
the collective behavior of the small- x gluons. The CGC approach uses classical color fields: 
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Fig. 2: Left: typical diagram for the production of high— pr particles, with large values of x being probed in the 
nuclear wave functions. Right: typical diagram for the production of bulk particles with pr ~ Q s , where multiple 
partonic interactions are crucial. This is true in heavy-ion collisions, and pp collisions at very high energies. 

The long-lived, large-x partons are represented by a strong color source p~ 1/gs which is static 
during the lifetime of the short-lived small-x gluons, whose dynamics is described by the color 
field A ^ 1/gs- The arbitrary separation between the field and the source is denoted xa- When 
probing the CGC with a dilute object carrying a weak color charge, the color field A is directly 
obtained from p via classical Yang-Mills equations: 



and it can be used to characterize the CGC wave function <& XA [A] . 

This wave function is a fundamental object of this picture, it is mainly a non-perturbative 
quantity, but the xa evolution can be computed perturbatively. Requiring that observables are 
independent of the choice of xa, a functional renormalization group equation can be derived. 
In the leading-logarithmic approximation which resums powers of as \ii(1/xa), the JIMWLK 
equation describes the evolution of \$ XA [A]\ 2 with xa- The evolution of the saturation scale with 
x is then obtained from this equation. 

Finally, the information contained in the wave function, on gluon number and gluon cor- 
relations, can be expressed in terms of n-point correlators, probed in scattering processes. These 
correlators consist of Wilson lines averaged with the CGC wave function, and resum powers of 
gsA ~ 1, i.e. scattering with an arbitrary number of gluons exchanged. In the CGC picture, both 
multiple scatterings and non-linear QCD evolution are taken into account. Note that in terms of 
occupation numbers, in the saturation regime one reaches 



Therefore, taking into account multiple scatterings in the collision is as important as the satura- 
tion of the wave function. A consistent calculation of MPI must include both. 

It was not obvious that the CGC picture ((U), which requires small values of xa, would be 
relevant at present energies. One of the most acclaimed successes came in the context of d+Au 
collisions at RHIC, where forward particle production pA — > hX allows to reach small values 
of xa with a dilute probe well understood in QCD [2]. The prediction that the yield of high-p^ 
particles at forward rapidities in pA collisions is suppressed compared to A pp collisions, and 
should decrease when increasing the rapidity, was confirmed. 
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Fig. 3: The charged-particle multiplicity in AA collisions at RHIC and the LHC. In both approaches a few parameters 
are fixed to reproduce RHIC data, such as the initial value of Q s . Then the small-a; evolution determines the multi- 
plicity at the LHC. The predictions are similar, around 1400 charged particles at mid rapidity for central collisions. 



3 Multiple partonic interactions in the Glasma 

The Glasma is the result of the collision of two CGCs. In a high-energy heavy-ion collision, 
each nuclear wave function is characterized by a strong color charge, and the field describing the 
dynamics of the small-x gluons is the solution of 

[D^, F^] = 5 +u pi + 8~ v pi . (5) 

The field after the collision is non-trivial [3]: it has a strong component {A^ ~ l/g s ), a compo- 
nent which is particle like (j4 m ~ 1), and components of any strength in between. To understand 
how this pre-equilibrium system thermalizes, one needs to understand how the Glasma field de- 
cays into particles. Right after the collision, the strong field component contains all modes. Then, 
as the field decays, modes with pr > 1/r are not part of the strong component anymore, and for 
those a particle description becomes more appropriate. After a time of order l/Q s , this picture 
breaks down, and it has been a formidable challenge to determine weather a fast thermalization 
can be achieved within this framework, due to instabilities [4]. 

A problem which can be more easily addressed is multiparticle production. The difficult 
task is to express the cross-section in terms of the Glasma field, and this is when MPI must be 
dealt with, as pictured in Fig.2. This has first been done at tree level, and from the one-loop 
calculation a factorization theorem could then be derived [5] (note an interesting possible appli- 
cation of the results to pp collisions: those first-principle calculations could inspire a model for 
the underlying event). Predictions for the total charged-particle multiplicity in AA collisions at 
the LHC are shown in Fig.3. Two approaches are compared: in the first, a simplified factoriza- 
tion (called kx factorization) is assumed but the energy evolution is accurately obtained from a 
next-to-leading evolution equation [6] ; in the second, the energy evolution is only parameterized 
but MPI are correctly dealt with by solving classical Yang-Mills equations [7]. While a full next- 
leading treatment of both multiple scatterings and small-x evolution is desirable, the numbers 
obtained are similar, which indicates that the uncertainties in both approaches are under control. 




Fig. 4: Left: production of high-energy partons in a hard process, which then lose energy propagating through the 
plasma. Some quantum fluctuations in their wave function are put on shell while interacting with the medium and be- 
come emitted radiation. Right: the resulting particle production in AA collisions is suppressed (Raa < 1) compared 
to independent nucleon-nucleon collisions. The suppression is large for light hadrons, and similar for heavy mesons 
(those data are displayed in the figure), which is difficult to accommodate in a weakly-coupled QCD description. 



4 The saturation scale in the QCD plasma 

Hard probes are believed to be understood well enough to provide clean measurements of the 
properties of the QGP formed in heavy-ion collisions. A large amount of work has been devoted 
to understand what happens to a quark (of high energy E, mass M and Lorentz factor 7 = E/M) 
as it propagates through a thermalized plasma [8]. MPI are a main ingredient of the perturbative 
QCD (pQCD) description of how a quark losses energy, until it thermalizes or exits the medium 
(see Fig.4). 

At lowest order with respect to a s , quantum fluctuations in a quark wave function consist 
of a single gluon, whose energy we denote lo and transverse momentum k± . The virtuality of that 
fluctuation is measured by the coherence time, or lifetime, of the gluon t c = u/k\. Short-lived 
fluctuations are highly virtual while longer-lived fluctuations are more easily put on shell when 
they interact. The probability of the fluctuation is a s N c , up to a kinematic factor which for heavy 
quarks suppresses fluctuations with u> > ^kj_. This means that when gluons are put on-shell, they 
are not radiated in a forward cone around a heavy quark. This suppression of the available phase 
space for radiation, the dead-cone effect, implies less energy loss for heavier quarks [9]. 

In pQCD, medium-induced gluon radiation is due to multiple scatterings of the virtual 
gluons. If, while undergoing multiple scattering, the virtual gluons pick up enough transverse 
momentum to be put on shell, they become emitted radiation. The accumulated transverse mo- 
mentum squared picked up by a gluon of coherence time t c is 

p± = i^ 2 j = qt c (6) 

where is the average transverse momentum squared picked up in each scattering, and I is the 
mean free path. These medium properties are involved through the ratio q = n 2 /l. 



Since only the fluctuations which pick up enough transverse momentum are freed (k± < 
p±), the limiting value can be obtained by equating k\ with p\ = q\o jk\: 

k± < {quj) 1/A = Qs(lo) . (7) 

The picture is that highly virtual fluctuations with k± > Q s do not have time to pick up enough 
p± to be freed, while the longer-lived ones with k± < Q s do. That transverse momentum Q s 
which controls which gluons are freed and which are not is called the saturation scale. With 
heavy quarks, one sees that due to the dead cone effect, the maximum energy a radiated gluon 
can have is u = ^k±_ = ~/Q s (and its coherence time is t c = j/Q 3 ). This allows to estimate the 
heavy-quark energy loss: 

- ^ oc a s N c ?%- = a s N c Q 2 s . (8) 
at j/Qs 

The saturation momentum in this formula is the one that corresponds to the fluctuation which 
dominates the energy loss: Q s = (qj) 1 ^ 3 - 

For a plasma of extend L <t c = 7 2 / 3 /g 1 / 3 , formula (O still holds but with Q 2 S = qL. 
These are the basic ingredients of more involved phenomenological calculations, but after com- 
parisons with data, it has remained unclear if this perturbative approach can describe the suppres- 
sion of high— p± particles. For instance, at RHIC temperatures, the value q ~ 1 — 3 GeV 2 /fm is 
more natural than the 5 — 10 GeV 2 /fm needed to describe the data on light hadron production. 
If one accepts to adjust q to this large value, then the D and B mesons are naturally predicted to 
be less suppressed than light hadrons, which is not the case (see Fig.4). 

While the present pQCD calculations should still be improved, and may be shown to work 
in the future, this motivated to think about strongly-coupled plasmas. The tools to address the 
strong-coupling dynamics in QCD are quite limited, however for the N = 4 Super- Yang-Mills 
(SYM) theory, the AdS/CFT correspondence is a powerful approach used in many studies. The 
findings for the strongly-coupled SYM plasma may provide insight for gauge theories in general, 
and some aspects may even be universal. One interesting result is that the total energy loss of 
hard probes goes as AE oc L 3 at strong coupling [10], instead of the L 2 law at weak coupling. 
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